Inequivalent routes across the Mott transition in V2O3 explored by X-ray absorption 
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The changes in the electronic structure of V2O3 across the metal-insulator transition induced by 
temperature, doping and pressure are identified using high resolution x-ray absorption spectroscopy 
at the V pre if -edge. Contrary to what has been taken for granted so far, the metallic phase reached 
under pressure is shown to differ from the one obtained by changing doping or temperature. Using 
a novel computational scheme, we relate this effect to the role and occupancy of the a\ g orbitals. 
This finding unveils the inequivalence of different routes across the Mott transition in V2O3. 

PACS numbers: 78.70.En; 71.27.+a; 71.20.Eh 
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Some materials present metal-insulator transitions 
(MIT) without any changes in crystal structure or long- 
range magnetic order. These phenomena, known as 
Mott-Hubbard transitions, constitute a fundamental sig- 
nature of strong electronic correlations. The physics 
emerging in the vicinity of these transitions is highly 
non-trivial and the properties of such materials depend 
on small changes in the electronic structure induced by 
external parameters [l], @] • Several features of the MIT 
have been successfully clarified by resorting to realistic 
many-body calculations 0. Yet, contrary to common as- 
sumptions, a growing number of experimental facts are 
revealing that this physical process is also strongly de- 
pendent on the route followed through the MIT. Here 
we show this to be the case for vanadium sesquioxide 
(V 2 3 ). The isostructural MIT in Cr-doped V 2 3 is 
considered as the textbook example of a Mott transition, 
which occurs between a paramagnetic insulator (PI) and 
a paramagnetic metallic (PM) phase by changing dop- 
ing level (x), temperature (T) or pressure (P) [4|. As 
such V2O3 has served as a test bed for many theoretical 
models [H4lH and a sustained experimental effort. 

Among the different experimental methods recently 
employed to study the electronic properties of the Mott 
transition in Cr-doped V2O3 ID" 15], X-ray absorption 
spectroscopy (XAS) has played a crucial role. It was the 
detailed investigation of the V L 2 ,3 absorption edges [l(| 
that demonstrated the necessity of abandoning the sim- 
ple one band, S = 1/2, model to obtain a realistic de- 
scription of the changes in the electronic structure at the 
phase transition. XAS can also be performed at the V 



K-edge in the hard x-ray range: in this case, the pre- 
edge will carry most of the physical information we are 
interested in, as it is predominantly due to Is — > 3d tran- 
sitions. The excitations in this spectral region are influ- 
enced by the core hole and should be considered to be of 
an excitonic nature. The advantages over the £2,3 edge 
studies are (i) a more straightforward interpretation (due 
to the simple symmetry of the s-core hole, the multiplet 
structure reveals a direct view on the d-states) and (ii) 
the possibility of applying external pressure as diamond 
anvil cells are compatible with hard X-rays. 

We used V K-edge XAS to explore extensively the MIT 
in V2O3 by changing temperature, doping and apply- 
ing an external pressure. The pre K-edges were ana- 
lyzed by a novel computational scheme combining local 
density approximation and dynamical mean field theory 
(LDA+DMFT) [l7| with configuration interaction (CI) 
full multiplet ligand field calculations to interpret subtle 
differences at the PM-PI transition. This allowed us to: 
(i) observe in detail the changes in the electronic exci- 
tations, providing also a direct estimate of the Hund's 
coupling J; (ii) analyze the physical properties of the PI 
and PM phase on both sides of the MIT, leading to the 
main result of our work: (iii) understand the difference 
between P, T or doping-induced transitions. This dif- 
ference is mainly related to the occupancy of the a\ g or- 
bitals, suggesting the existence of a new "pressure" path- 
way between PI and PM in the phase diagram. 

The experiments were performed on the inelastic x-ray 
scattering (IXS) beamlines ID-26 and BM-30 at the Eu- 
ropean Synchrotron Radiation Facility (ESRF). We used 
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high quality samples of (Vi-xCr^C^ with various dop- 
ing in the PM (x = 0) and PI phases (x = 0.011, and 
0.028) at ambient conditions. The MIT was also crossed 
for the 0.011 doping by changing temperature and for 
the 0.028 doping by pressure. To obtain the best reso- 
lution, the XAS spectra were acquired in the so-called 
partial fluorescence yield (PFY) mode [HI, monitoring 
the intensity of the V-Ka (2p — s> Is) line as the incident 
energy is swept across the absorption edge. Compared to 
standard XAS, the PFY mode provides better resolved 
absorption spectra as these are partly free from core-hole 
lifetime broadening effect. The gain is particularly ap- 
preciable at the transition metal K pre-edges as the Is 
core-hole is extremely short lived with respect to the 2p 
core-hole and the dipolar tail is strongly reduced. The 
V-Ka line was measured with the help of a IXS spec- 
trometer equipped with a Ge(331) spherically bent ana- 
lyzer. For the temperature dependence, the crystals were 
mounted in a cryostat installed on the spectrometer. On 
the other hand, measuring the PFY mode turned out to 
be difficult under pressure because of the weak fluores- 
cence signal in the pressure cell and we opted for standard 
XAS measurements in transmission mode. To maximize 
the throughput, a powder sample (x — 0.028) was loaded 
in a diamond anvil cell equipped with composite anvils 
made of a perforated diamond capped with a 500-/im 
thin diamond anvil. Pressure was measured in-situ by 
standard ruby fluorescence technique. We used silicon 
oil as a pressure transmitting medium. 

The T-dependent absorption spectra are displayed in 
Fig. [Ha) for both PM (200 K) and PI (300 K) phases 
for the x = 0.011, powder sample. The spectra have 
been normalized to an edge jump of unity. We will focus 
on the pre-edge region, where information about the V d- 
states can be easier extracted. It can be decomposed into 
three spectral features (A,B,C) which all vary in intensity 
as the system is driven through the MIT whereas only 
C is considerably broadened (Fig. d{a), inset). Notice 
that no feature is observed below peak A contrary to the 
early results of Ref. 11911 but in agreement with the more 
recent data of Rcf . 20] . Within a simplified atomic like 
picture, one could directly relate the intensity of features 
A,B and C to the unoccupied states. Starting from a V- 
t^g, S — 1 configuration, one can either add an electron 
to the tig subshell (split into one ai g and two e£ states 
under trigonal distortion of the V sites (HJ as shown in 
Fig. [3Jd)) yielding peaks A and B, or add an electron 
to the sub-shell which gives rise to the broader peak 
C. In this picture, one could associate peak A with a 
spin quartet (S — 3/2) and peak B with a spin doublet 
(S = 1/2), split by the Hund's rule exchange. 

This picture is however oversimplified as the V d elec- 
trons are strongly correlated and the spectra are still 
largely influenced by the Is core hole potential. Keeping 
that in mind, we have simulated the pre-edge by combin- 
ing CI with LDA+DMFT calculations for which the one 
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FIG. 1: (Color online) (a) V if-edge x-ray absorption spec- 
tra in (Vi-^Cr^Os, x = 0.011 powder sample measured 
as a function of temperature (T) in the PM (200 K, solid 
line) and PI (300 K, dotted line) phases by PFY XAS; below, 
PM— PI spectral difference; (inset) pre-edge region, (c) Cal- 
culated isotropic CI V K-e&ge PM and PI spectra of V2O3. 
The spectra are compared to the T-dependent experimental 
data shown in the same energy window (b). (d) From LDA 
(starting point for our LDA+DMFT and CI calculations) we 
obtain that the cubic part of the ligand field splits the V d- 
levels in ti g and e g by lQDq ~2 eV, and that the small trigonal 
distortion further separates the t^ g into a\ g and e g by Atrig 
-0.3 eV. 



particle part (LDA) input corresponds to the level dia- 
gram in Fig. []Jd) . We concentrate our analysis to peaks 
A and B, since peak C relates mainly to the unoccupied 
e a g orbitals. These hybridize much stronger with the oxy- 
gen ligands and thus lack direct information on the Mott 
transition; peak C may also be related to non-local exci- 
tations (not included here) [22} which sensitively depend 
on the metal-ligand distance. Let us also note that the V 
sites in V2O3 are non centrosymmetric which leads to an 
on-site mixing of V-3ri and V-4p-orbitals and interference 
between dipole and quadrupole transitions 17] which has 
been included in our scheme. 

Our CI calculations confirm that for the ground state 
the occupancy ratio between the (eg,ai g ) and (e![ ,e" ) 
states is smaller in the PI than in the PM phase [81. 1 161] : 
The isotropic Cl-based calculated XAS spectra in the 
pre-edge region reported in Fig. [He) agree well with the 
experimental data for both the energy splitting of fea- 
tures A and B and the ratio of their spectral weight (SW) 
which increases in the PM phase. 
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Considerable insight can be gained by comparing CI 
and LDA+DMFT calculations. Our LDA+DMFT cal- 
culations, performed using N t,l -order muffin-tin orbital 
(NMTO) downfolded Hamiltonian for the 1.1% Cr-doped 
V2O3 and Hirsch-Fye Quantum Monte Carlo as impurity 
solver, confirm the above mentioned tendency (we obtain 
a mixing of 50:50 and 35:65 for the (e£ ,a\ g ):(e^ ,ep oc- 
cupation in the PM and PI phases). For LDA+DMFT 
we have chosen the same interaction parameter U = 4.2 
eV for the PI phase as in Ref. [§] and followed their sug- 
gestion to slightly decrease its value in the PM phase (we 
assume U = 4.0 eV)j23[. Remarkably the simple struc- 
ture of the core hole potential in the K-edge spectroscopy 
allows us to associate the pre-edge spectrum with the 
k— integrated spectral function above the Fermi energy 
calculated by LDA+DMFT. The electron-addition part 
of the spectral function shows three main features in PM 
phase: a coherent excitation at the Fermi level and a 
double peak associated to the incoherent electronic ex- 
citations i.e. the upper Hubbard band (UHB), similarly 
to the undoped compound. In the PI phase obviously, 
only the latter survives. Comparison with the experi- 
mental spectra clearly shows that the pre-edge features 
have to be related to the "incoherent" part of the spectral 
function only. The physical reason is that the core hole 
potential localizes the electrons destroying the coher- 
ent quasiparticle excitations and making the XAS spec- 
trum atomic-like. All the "incoherent" LDA+DMFT, 
CI, and experimental spectra shown in Fig. [2] agree in 
many aspects, especially as for the splitting of the first 
two peaks by «2.0 eV («1.8 eV in experiment) which 
originates in LDA+DMFT from the Hund's exchange J 
in the Kanamori Hamiltonian. This further validates 
the choice of J = 0.7 eV used in our calculations con- 
trasting with larger values assumed in previous stud- 
ies [1, [13, HH, and also clarifies the mismatch between 
XAS and LDA+DMFT spectra reported in the undoped 
V2O3 compound 0]. Moreover, the ratio between A and 
B peak displays the same trend in the PM-PI transition 
as the CI (or experimental) data. The quantitative dif- 
ference between the two calculations is attributed to the 
lack of matrix elements in LDA+DMFT. 

The intensity ratio of the first two incoherent excita- 
tions peaks A and B (associated to the quartet and dou- 
blet states in the oversimplified picture) thus appears as 
the key spectral parameter to understand the differences 
between PM and PI. Even in a powder sample, this ratio 
is still sensitive to the a\ g orbital occupation of the initial 
state. Indeed, due to the trigonal distortion a consider- 
able spectral weight transfer from the peak B to higher 
energies (corresponding to final states with two a\ g elec- 
trons in the limit of large A tr i g ) can take place for the 
(e^, ai g ) but not for the [e K g ,e K g ) initial state. Therefore, 
the K pre-edge XAS can serve as a direct probe of the 
aig orbital occupation in the ground-state. As a rule of 
thumb, the larger the ratio between the SW of A and B, 
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FIG. 2: (Color online) Incoherent LDA+DMFT and CI cal- 
culations in the pre-edge region; Ef is the Fermi energy. Note 
the similarity in the main spectral features when crossing 
the MIT. Also shown are the different contributions of the 
CI spectrum labeled accordingly to their initial state: the 
contribution of the (e^,a\ g ) — > (eZ ,el ,ai 3 ) transitions to the 
peak A(B) is approximately 60%(55%) in the PM phase and 
20%(15%) in the PI phase. 

the larger the a\ g orbital occupation. 

Hard X-ray absorption also provides a unique way to 
explore the "pressure" pathway across the MIT, from 
which relevant information can be extracted by applying 
our interpretation scheme. Fig.|3][a) shows the XAS pow- 
der spectra of the P-induced MIT with the corresponding 
spectra for the T- and doping-driven transition (cf. the 
loci in the phase diagram, Fig. [3]). To ease the com- 
parison, the pressure dataset (P) obtained in transmis- 
sion mode has been deconvolved from the Is Lorentzian 
lifetime broadening (1 eV FWHM) using the GNXAS 
code [2^ to match the improved resolution of datasets 
(x) and (T) measured by PFY-XAS. Together with the 
PM— PI differences in Fig. [2Kb) , Fig. 131(a) clearly evi- 
dences that, besides a rigid shift of the first two peaks of 
~+0.13 eV and contrary to the doping- or T-driven tran- 
sition, very small changes in spectral shapes and weights 
are observed in the P-driven MIT. In the light of the 
arguments discussed above, our finding proves that the 
metallic state reached by applying pressure is character- 
ized by a much lower occupation of the a\ g orbitals com- 
pared to the metallic state reached just by changing tem- 
perature or doping. Importantly, the spectra measured 
through the doping induced MIT are identical within 
the experimental uncertainty to those measured through 
the T-driven transition. The temperature-doping equiv- 
alence is confirmed by our photoemission data displayed 
in Fig. |3tc,d) [l5| and is consistent with the very similar 
lattice parameter changes across the transition [25(] . The 
x and T equivalence is also borne out by the observa- 
tion from XAS at the L-edges in doped V2O3 [13] that 
the a\ g occupation within both the PM or PI phases is 
mostly independent of the doping level. Hence, the local 
incoherent excitations probed by XAS at the V L edg e 
or K pre-edge are not directly affected by disorder [261 ] . 
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FIG. 3: (Color online) (a) V-K XAS spectra for powder sam- 
ples of (Vi-xCr^Oa as a function of pressure (x) [x — 0.028; 
5 and 7 kbar (lower curves); 5 and 11 kbar (upper curves)] 
(P), temperature (A) [x = 0.011; 200,300 K] (T), and dop- 
ing (o) [x = 0, 0.011] (x) (cf. points in the phase diagram; 
the pressure scale refers to the x = 0.028 doping). The spec- 
tral differences PM— PI (x3) shown for the three datasets (b) 
demonstrate the nonequivalence between P and temperature- 
doping. The x-T equivalence is confirmed by the photoemis- 
sion spectra (hu = 9 eV) (c) and their spectral differences (d) 
(from Ref. [3]). 



Our finding clearly shows the limits of the common as- 
sumption that temperature/doping- and pressure-driven 
MIT in V2O3 can be equivalently described within the 
same phase diagram [lajiliS]- Indeed, the two different 
PM electronic structures that we observed reflect differ- 
ent mechanisms driving the MIT along different path- 
ways. In the doping-driven MIT, the metallic phase 
is characterized by an increased occupation of the a\ g 
electrons indicating a reduced "effective crystal-ficld- 
splitting" as the main driving mechanism towards metal- 
licity [<3, 9], related to the jump of the lattice parameter 
c/a (1.4%) at the MIT [25|. In contrast, when pressure is 
applied the a\ g occupation remains basically unchanged, 
so that this metallic phase seems to originate rather from 
an increased bandwidth, without any relevant changes of 
the orbital splitting. The smaller c/a jump observed un- 
der pressure (0.7 %) corroborates our analysis. 

In conclusion, doping, temperature and pressure are 
shown to act differently on the interplay between electron 
correlations and crystal field, so that states previously 



considered to be equivalent metals are actually differ- 
ent. We believe this finding may apply to many other 
strongly correlated systems presenting metal-insulator 
transitions, with pressure opening inequivalent pathways 
through their phase diagrams. 

We thank E. Tosatti for stimulating discussions and ac- 
knowledge a BQR of the Universite Paris-Sud, the "Tri- 
angle de la Physique", the Austrian Science Fund (FWF, 
science college W004), and the EU research network RP7 
MONAMI for financial support. 



[1] M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 
70, 1039 (1998). 

E. Dagotto, Science 309, 257 (2005). 
A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

D. B. McWhan, A. Menth, J. P. Remeika, W. F. 
Brinkman, and T. M. Rice, Phys. Rev. B 7, 1920 (1973). 

C. Castellani, C. R. Natoli, and J. Ranninger, Phys. Rev. 
B 18, 4945 (1978). 

K. Held, I. Nekrasov, N. Bliimer, V. Anisimov, and 

D. Vollhardt, Int. J. Mod. Phys. B 15, 2611 (2001). 
I. S. Elfimov, N. A. Skorikov, V. I. Anisimov, and G. A. 
Sawatzky, Phys. Rev. Lett. 88, 015504 (2002). 

G. Keller, K. Held, V. Eyert, D. Vollhardt, and V. I. 
Anisimov, Phys. Rev. B 70, 205116 (2004). 
A. I. Poteryaev, J. M. Tomczak, S. Biermann, 
A. Georges, A. I. Lichtenstein, A. N. Rubtsov, T. Saha- 
Dasgupta, and O. K. Andersen, Phys. Rev. B 76, 085127 
(2007). 

M. S. Laad, L. Craco, and E. Muller-Hartmann, Phys. 
Rev. B 73, 045109 (pages 15) (2006). 
A. Tanaka, J. Phys. Soc. Jpn. 71, 1091 (2002). 
P. Limelette, A. Georges, D. Jerome, P. Wzietek, P. Met- 
calf, and J. M. Honig, Science 302, 89 (2003). 
S.-K. Mo, J. D. Denlinger, H.-D. Kim, J.-H. Park, J. W. 
Allen, A. Sekiyama, A. Yamasaki, K. Kadono, S. Suga, 
Y. Saitoh, et al., Phys. Rev. Lett. 90, 186403 (2003). 
S.-K. Mo, H.-D. Kim, J. D. Denlinger, J. W. Allen, J.- 

H. Park, A. Sekiyama, A. Yamasaki, S. Suga, Y. Saitoh, 
T. Muro, et al., Phys. Rev. B 74, 165101 (2006). 

F. Rodolakis, B. Mansart, E. Papalazarou, S. Gorovikov, 
P. Vilmercati, L. Petaccia, A. Goldoni, J. -P. Rueff, 
S. Lupi, P. Metcalf, and M. Marsi, Phys. Rev. Lett. 102, 
066805 (2009). 

J.-H. Park, L. H. Tjeng, A. Tanaka, J. W. Allen, C. T. 
Chen, P. Metcalf, J. M. Honig, F. M. F. de Groot, and 

G. A. Sawatzky, Phys. Rev. B 61, 11506 (2000). 
K. Held, Adv. in Phys. 56, 829 (2007). 
F. M. F. de Groot, Chem. Rev. 101, 1779 (2001). 
A. Bianconi and C. R. Natoli, Solid State Commun. 27, 
1177 (1978). 

J. Goulon, A. Rogalev, C. Goulon-Ginet, G. Benayoun, 
L. Paolasini, C. Brouder, C. Malgrange, and P. A. Met- 
calf, Phys. Rev. Lett. 85, 4385 (2000). 
T. Saha-Dasgupta, O. K. Andersen, J. Nuss, 
A. I. Poteryaev, A. Georges, and A. I. Lichtenstein, 
arXiv:0907.284Tl (preprint 2009). 

C. Gougoussis, M. Calandra, A. Seitsonen, C. Brouder, 



[2 
[3 

[4; 

[5 

[6 
[7; 
[8 
[9 

[10: 

[11 
[12 

[13 
[14 
[15 

[16 



[17 
[18 
[19 

[20 



[21 



[22 



5 



A. Shukla, and F. Mauri, Phys. Rev. B 79, 045118 (2009). 

[23] Constrained LDA/LDA+U estimates of J (0.93eV in 
[28(| ) need to be reduced to J ~ 0.7eV to reproduce 
the constrained LDA energy splitting with the DMFT 
Hamiltonian [vf\ . Moreover, the value of U adopted for 
LDA+U (2.8eV) has to be assumed smaller than that for 
LDA+DMFT to compensate deficiencies of LDA+U [2<|. 

[24] A. Filipponi, J. Phys. B: At. Mol. Opt. Phys. 33, 2835 
(2000). 

[25] D. B. McWhan and J. P. Remeika, Phys. Rev. B 2, 3734 
(1970). 

[26] A. I. Frenkel, D. M. Pease, J. I. Budnick, P. Metcalf, E. A. 



Stern, P. Shanthakumar, and T. Huang, Phys. Rev. Lett. 

97, 195502 (2006). 
[27] An early version of the phase diagram (Fig. 15 in Ref . I25T) 

was drawn with a third pressure axis: due to the idea of 

P-x equivalence, this was later abandoned. 
[28] I. Solovyev, N. Hamada, and K. Terakura, Phys. Rev. B 

53, 7158 (1996). 
[29] G. Sangiovanni, A. Toschi, E. Koch, K. Held, M. Capone, 

C. Castellani, O. Gunnarsson, S.-K. Mo, J. W. Allen, H.- 

D. Kim, et al., Phys. Rev. B 73, 205121 (2006). 



